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6-SEP=1 
-TITLE MTHSDTANH 


2 
2 fa ; Floating Point Hyperbolic Tangent rout ae: EP=1 138 
+ 
IDENT /1-011/ ; F 


—BARAAAAAARRAAAALALALAALALALAL EAA ASE EERE EAA R SESE REE ER REAR REE RR SEES SE SY 


COPYRIGHT (c) 1978, 1980, 1982, 1984 B 
DIGITAL EQUIPMENT CORPORATION, MAYNARD. MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


® 
® 
® 
* 
® 
® 
THIS SOFTWARE 1S FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED 
ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE * 
INCLUSION OF THE ABOVE COPYRIGHT NOTICE. THIS SO FTWARE OR ANY. OTHER 
COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY * 
OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY * 
TRANSFERRED. * 
oe 
® 
® 
® 
ics 
® 
® 
wt 
® 
® 


THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
preeokat ihe NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 


COOCOOCOCCOOCOOOOCOoOoOoOooooo 
COOQooooooQooQooQooooooooo 


SOOoQooooooooooo 
COOQCoooooooooooooooooo 
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19 ; 
11 ; 
i 3 
2 8 
14 ; 
12 3 
16 ; 
i 3 
18 ; 
19 ; 

0000 3* DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 

8000 : SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 

0000 ; 

0000 3 REARAEAAAEAAAATATAAETAAAKAAEAAAAAAARAAAAAAAAAAAAAARAAAAEAEARAEARAARAAAAAARAERAEARERAEEAEE 

0000 3 

0000 $ 

0000 3; FACILITY: MATH LIBRARY 

0000 gee 

+44 3 ABSTRACT: 

0000 3; MTHSDTANH is a function which returns the floating point hyperbolic tangent 
000 3; of its single precision floating point argument. The call is standard 
000 3 call-by-reference. 

0000 3 

0000 ;-- 

0000 3 

0000 39 ; VERSION: 01 

a4 40 ; 

000 41 ; eraerests 

0000 42 ; AUTHOR: : 

444 2 : Peter Yuo, 29-Jun-77: Version 01 

6009 45 ; MODIFIED BY: 

000 46; 

4 47 ; 
000 48 ; 
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1-011 liStORY e pecet 
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5 
4 
BY 
00 $ 
$608 
000 60 
$00 
$00 
0000 64 
0000 «6S 
0000 68 
0000 «6? 
0000 68 
000 69 
0000 70 
0000 «71 
9000 72 
0000 «7 
0000 «7% 
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urrent Edit History 6-SEP-1984 MTHRTL SRE SATHDTANH. MAR: 1 
-SBTTL HISTORY ; Deta:icd Current Edit History 


ALGORITHMIC DIFFERENCES FROM FP-11/C ROUTINE: none 
Edit History for Version 01 of MTHSDTANH 


=3 Call getty a directly instead of using POLY. 
=4 Add .EXTRNs. TNH 13-June-7 
5 - Use W sitent on externals. TNH 13-June-78 
More accurrate. TNH 13-June-78 
1 06 - Update pares pusher and copyright Eres ices jBs 16-NOV-78 
1-007 = Add ‘'"' to the PSECT directive. -7 
1-008 - Use MTHSDEXP_R6. SBL 37e§ te t-197 
1-009 = Change constant 16.0 to to correct inaccuracy. 
value of X above which #¢ is the best machine eros 
to DTANH(X) is about 20.10. The neat higher number that can 
be represented as a short literal is et § JAW 19-Sep-80 
1-010 = Use general mode aera: = L_30-Nov-198 
1-011 - Natge constant 2*-1 *-28 to correct inaccuracy. vor 
es of iX! between vs -i2 and 2*-28 DTANH was only accurate 
a ig te hy since the soeunes ton that ag pee g in that 
of values is false. DTANH(X)= for ix! {ene -2 
DTA H(XD= DSINH(X) /DCOSHEX) for 2°-28<iXi<=. - All appropriate 
references to 2°-14 have been changed to 33238" JCw Jan-1983 
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1-011 BEC RRET TONS obec Loree ee Pe eee MCCUE GxSE patos $1 :53057 PATANTS caecariorana.man:1 28 cd) 1- 
re -SBTTL DECLARATIONS : Declarative Part of Module 
hi 3 INCLUDE FILES: 
1 : EXTERNAL SYMBOLS: MTHSJACKET_HDLR 
° é -DSABL GBL 3 Force .EXTRN on all symbols 
4 -EXTRN MTHSDCOSH ; DCOSH 
0 5 -EXTRN MTHSDSINH ; DSINH 
0 6 “EXTRN MTHSDEXP_R6 : EXP 
8 Ms 3; EQUATED SYMBOLS: 
0004080 000 9 $0_1.0 = *F1.0 3: 1.0 
0004280 000 0 $D722.0 = *F22.0 3 22. 
00000004 000 91 value =4 3 value.rd.r 
000 3 
000 93 ; 
44 94 ; MACROS: none 
00 7) 3 
444 38 3 PSECT DECLARATIONS: 
00000000 98 -PSECT _MTHSCODE PIC,SHR,LONG, EXE ,NOWRT 
44 190 3 program section for math routines 
0000 +2 3; OWN STORAGE: none 
000 1 ¢ ; 
0000 103; 
0000 104 ; CONSTANTS: 
0000 105; 
0000 106 
0000 107 D_0.25: 
0000 0000 0000 3F80 0000 198 «WORD “xX3F80, 0, 0, 0 : 0.25 
0008 109 D_2_POWER M28: 
0000 0000 0000 3280 44 19 -@ORD *x3280, 0, 0, 0 3 2ee-28 
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-SBTTL MTHSDTANH = Standard DOUBLE Precision Floating DTANH 


“FUNCTIONAL DESCRIPTION: 
DTANH = double precision floating ooint function 
DTANH(X) is computed as: 

If ik! =< 2e*-28, then DTANH(X) = X. 


If 2e*-28 < IX! =< 0.25, then DTANHiX) = DSINH(X)/DCOSH(X). 
If 6.25 < 1X! < 22.0; then DTANH(X) = (DEXP(2*x) = 1) / (DEXP(2*x) + 1) 
If 22.0 =< IX}, then DTANH(X) = sign(x) * 1 


CALLING SEQUENCE: 
DTANH.wd.v = MTHSDTANH(x.rd.r) 
INPUT PARAMETERS: 


LONG = 4 :; define songeers multiplier 
x = 1 * LONG 3; Contents of x is the argument 


IMPLICIT INPUTS: none 
OUTPUT PARAMETERS: 
VALUE: double precision floating hyperbolic tangent of the argument 
IMPLICIT OUTPUTS: none 
COMPLETION CODES: none 
SIDE EFFECTS: none 


NOTE: This procedure disables floating point underflow, enables integer 
overflow. 


ENTRY MTHSDTANH, “M<IV, R2, R3, R4, RS, ROD 

3 standard call-by-reference entry 

3; disable DV (and FU), enable IV 
MTHSFLAG_JACKET 3; flag that this is a jacket procedure in 


MOVAB G*MTHSSJACKET_HND, (FP) 
3 set handler address to jacket 


; handler 
; case of an error in routine 
:; If an error, convert signal to user PC 
; and resigna 
mova Bugging MOE ARAL Sy = Batwa 
CMPD = RO, Do POMER_M28 > compare !X! with 22-28 
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00000000' 
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00000000 ' GF g 
00000000°GF 6 
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andard DOUBLE Precision 
BLEQ OUT_X 


2ee-28 < IX! 


CMPD 
BGEQ 


> 2ee-28 < IX! < 22.0 


CMPD 
BLEQ 


: 0.25 < ix! < 22.0 


2 DeeeIR6 < IX! =< 0.25 


LEQ_T0_0.25: 
CALLG 


MOVD 

CALLG 

DIVD 

RET 
GEQ_T0_22.0: 

MOVD $*#SD_1 

TSTD apgtues 

BGEQ 10$ 

MNEGD RO, RO 
10$: RET 


Geé_to.2270 


RO, 0.0.25 
LEG_10_0.25 


@value(AP), a@value(AP), 


ADDD3 

JSB G*MTHSDEXP_R6 
ADDD3 $*#SD_1.0 
SUBD  $*#SD-1.0 
DIVD R2, RO 

RET 


(ap) G“MTHSDCOSH 
(AP), G*MTHSDSINH 
R2, RO 


@value(AP), RO 


AX/VMS Macro V04-00 
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3; branch if {Xi =< 2*#-28 


Sete te te te te 
Oo 


compar 
branch 


compar 


*t int Se 35° 


e ix! with 0.33 


branch if ix: =< 0. 


return w 


RO/R1 
test t 
branch 
RO/R 
return 


RO/R1 
return 


wer 
‘+ 
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esult in RO/R1 


DCOSH(X) 
DCOSH(X) 
DSINH(X) 
DSINH(X) / DCOSH(X) 
ith result in RO 


= 1.0 
he sign of X 
if X >= 0 


= 
with result in RO 


= DTANH(X) = X 
with result in RO/R1 
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1 
(2X) = 1) / (DEXP(2*xX) + 1) 
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MTHSDTANH 3 Floating Point Hyperbolic Tan A rout 16-SEP-1984 01: AX/VMS M 
Symbol table . ‘nt Bent rout tee Tempe $1:53:59 YatANte ee 
8 1 
b-9 28 POW 8 M28 0000008 1 
GEA. -10. 250 FR 1 
2 aA R 1 
= 
Lyte sacKeT HND Seseeeae x 1 
MTHSDCOSH reeeerese x 0 
MTHSDEXP_R6 eteenere =X 3 
MTHSOSINA rereeeee x 0 
MTHSDTANH 00000010 RG 01 
OUT Xx 000 eit, R 01 
$0_T.0 = 4494 4 58 
$D722.0 = 00 
VACUE = 90000 Be 
teecee meme eren aac + 
: ; Psect synopsis H 
PSECT name Allocation PSECT No. Attributes 
- ABS. 00000000 ( 9.) 0.) NOPIC USR CON ABS” LCL 
_MTHSCODE 00000070 ¢ 11 hel PIC USR CON REL LCL SHR_- EXE 
ened Bea 
Phase Page faults CPU Time Elapsed Time 
Initialization 9 0: 8:00.08 00:00:00.70 
zomnene processing 124 0:00: 8: 4 8:08: 5.34 
3 00:00:00.78 0:00: et 
Symbol table sort 9 Bp: : 2} Hf 01 
Pass 2 5 :00:00. 0:00:01.77 
Symbol table output 00:00:00.0 0:00:00.06 
Psect synopsis output 1 208: x 0: 208-8 
Cross-reference output 9 :00: $. 0 0:00:00.00 
Assembler run totals 29 :00:02.26 0:00:10.21 


ioe working rs asage was Pao ages. 
9 bytes ‘ pages) of virtua memory were used to buffer the intermediate code 
dl were 10 pages of symbol table space allocated to hold 15 non-local and 1 local symbols. 
84 source (anes were read in Pass 1, producing 11 object records in Pass 2. 
page uf virtual memory was used to define 1 macro. 


Geoceoeerccesoscesosesooose$ 


Macro Library name Macros defined 
~8255$DUA28: CSYSLIBISTARLET.MLB; 2 0 
0 GETS were required to define 0 macros. 


There were no errors, warnings or information messages. 
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se NOEXE NORD NOWRT NOVEC BYTE 
RD NOWRT NOVEC LONG 
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r 
SEP= L.SRCJMTHDTANH.MAR; 1 


MACRO/ENABLE=SUPPRESSION/DISABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHDTANH/OBJ=OBJ$:MTHDTANH MSRCS$:MTHJACKET/UPDATE=(ENHS$ :MTHJACKET) +MSRC 
, 


NT CORPORATION 
D PROPRIETARY 
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